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ABSTRACT 



We represent the partition function of the Generalized Kontsevich Model 
(GKM) in the form of a Toda lattice r-function and discuss various implica- 
tions of non-vanishing "negative"- and "zero"-time variables: the appear to 
modify the original GKM action by negative-power and logarithmic contribu- 
tions respectively. It is shown that so deformed r-function satisfies the same 
string equation as the original one. In the case of quadratic potential GKM 
turns out to describe forced Toda chain hierarchy and, thus, corresponds to a 
discrete matrix model, with the role of the matrix size played by the zero-time 
(at integer positive points). This relation allows one to discuss the double- 
scaling continuum limit entirely in terms of GKM, i.e. essentially in terms of 
finite-fold integrals. 
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1 Introduction 



In [1] a new matrix model was introduced which interpolates between all the non-perturba- 

Q(q-q') 2 

tive partition functions of Virasoro (q, g')-minimal string models with c = 1 . 

qq' 

The partition function of GKM depends on two distinct sets of time- variables: one entering 

s X p 

the "potential" V{X) = — — , and the other one is connected by a kind of Miwa 

p>\ P 

relation with auxiliary matrix M: 

t p = ^{TrM-v+( P -l)s p }. (1.1) 

For the particular choice of parameters (potential V(X) = const ■ X q+1 ) it coincides 
with the description of (q, g')-series£]. Interpolation as given by GKM preserves both nice 
properties of the known non-perturbative partition functions: if considered as functions 
of t-variables these are usually r-functions of integrable KP hierarchy and usually satisfy 
"generalized string equation" L^[\ = 0. These facts make GKM an appealing object to 
study in the context of non-perturbative and unified string theory. 

This paper describes some new enlightening results about GKM. The question to be 
addressed is what is the meaning of GKM from the point of view of the Toda hierarchy. 
The point is that the KP hierarchy can be always embedded into the Toda lattice hierarchy, 
but the latter one has additional variables: "negative times" and "zero-time" , which are 
frozen in KP case. The natural thing to ask is whether GKM can be further generalized to 
include the dependence on these additional parameters so that non-perturbative partition 
function becomes r-function of the Toda lattice hierarchy. 

Remarkably it appears possible not only to introduce the extra "zero-time" variable 
into GKM, but this is also a natural way to describe discrete matrix models, thus giving 
a chance to unify the entire theory of matrix models under a common roof of GKM. 

The natural form in which the non-perturbative partition functions^ arise, when de- 

1 E.g. for q = 2 - with Witten's topological gravity [2], which has been represented in the form of a 

matrix model by Kontsevich [3], 

2 We use the term "non-perturbative partition function" for the generating functional of all the exact 

[i.e. non-perturbative or appropriately summed over all orders of perturbation theory) correlation func- 
tions in string models. This object may be considered as a vacuum amplitude in the theory with the 
action, to which all the possible vertex operators (corresponding both to naive observables and to handle- 
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duced from the matrix models, is determinant of N x N matrix, probably with some extra 
simple (normalization) factors in front of it. Integrable r-functions are also representable 
through determinant formulae, but somewhat different expressions are adequately de- 
scribing different hierarchies. This representation of KP r-function arises naturally in 
Miwa parameterization of time-variables t p , i.e. essentially in terms of eigenvalues {rrii} 
of the N x N matrix M, and looks like 

{v} _ det &(mj) 
Tkp [tl ~ A(m) ' {l - 2} 

where {SAm) = m t ~ 1 (l + Oi — ))} is a Segal- Wilson basis, describing some point in the 

m 

Sato's Grassmannian (see all the details and references in [1]). On the other hand, the 
natural determinant formula for the Toda hierarchy is (see Appendix A, eq.(A.25)) 

r n [t] = det H i+nJ+n [t], (1.3) 
where satisfy the following equations: 

dHij/dtp = Hij- P for "positive times" t p , 

(1-4) 

dHij/dt_ p = Hi- P j for "negative times" t_ p , 
and n is integer- valued "zero-time" . 

In what follows we need more statements from the Toda theory. The forced Toda- 
lattice hierarchy arises if 



T -n[t] = $n,o for all n > 0. (1.5) 

The reduction of general Toda lattice hierarchy to the Toda chain one occurs, for example, 
when Hij depends only upon the difference i—j: Hij = Hi-^. (More details and references 
gluing operators) are added with arbitrary coefficients (which are nothing but the time- variables). Such 
quantity, though a priori constructed for one particular string model, a posteriori is naturally describing 
a family of models as large as the freedom, allowed in the deformations of the action. For a sufficiently 
rich family of deformations the quantity, describing the entire string theory (i.e. unification of all the 

string models), may be derived starting from any particular string model. 

3 The relations of this type, having a sense under the determinant, should be understood up to lower 

or upper triangle transformations admissible in the determinant. 
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concerning Toda hierarchies can be found in [4], see also (2.20) below). 

The partition function of GKM has been already proved in [1] to be a KP r-function 
and can be explicitly represented in the form of (1.2). Our first goal below will be to 
bring it to the form of (1.3). Usually KP hierarchy can be embedded into Toda lattice 
hierarchy at fixed values of zero- and negative-times^ 

T n [t. p ,t p } = 4f- p} [t p } . (1.6) 

In other words, the point of Grassmannian associated with the KP r-function at the r.h.s. 
depends on t_ p . This is also clear from the free-field representation of r-functions (see 
Appendix A): a Toda-lattice r-function, 

ri 9} [t- p ;t p ] ~< n\e~^ tpJp g e^'-^-p \ n > (1.7) 
can be considered CIS db KP r-function, 

T$*-> } [t p ] ~< n\e-^ Jp g KP \n > (1.8) 

with t_ p -dependent gxp = g e~^ t - pJ ~ p . In GKM the point g is specified by potential 
V(X) (this relation is conventionally encoded in the form of the "string equation"). We 
shall see that n- and /^-dependencies, as defined by eqs.(1.6) and (1.8), can be imitated 
by additional logarithmic and negative-power terms in the action of GKM. This is a sort 
of deformation of GKM in the sense that the original string equation [1] still remains 
to be valid. A natural question is what are the restrictions on GKM, which imply the 
occurrence of the forced Toda hierarchy. Essentially these requirements are for all the 
functions <fii(m) to be polynomial in m (up to trivial factor, see, for example, eq.(3.24)), 
and this is the case if V(X) = X 2 /2. The GKM as defined in [1] (i.e. with n = t_ p = 0) 
We systematically use the notation [. . .] to denote the dependence of a family of arguments like 
eigenvalues of M or all the time- variables (negative or positive). However, when the same objects are 
considered as parameters we put them into braces: {■••}• For example, negative times t- p are arguments of 
Toda lattice r-function, but they are parameters of KP r-function (they define the point in Grassmannian, 
but are not the time variables which are conventionally considered as arguments of tkp)- The main 
purpose of the study matrix models is certainly to get rid of this delicate difference, i.e. to treat all the 
parameters (which label different string models) as arguments. 
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is trivial for such potential, but as we shall see it becomes more interesting at n ^ and 

The second big problem addressed in this paper is the representation of discrete matrix 
models in the form of GKM. One should suggest that such representation exists if we 
indeed want GKM to be a kind of a universal matrix model (encoding all the information 
about the KP hierarchy and Grassmannian, which is relevant for string theory). This 
problem is also related to the Toda theory, since partition functions of discrete matrix 
models are known to be r-functions of forced Toda-lattice hierarchy (or even its Toda- 
chain reduction in the one-matrix case) [6]. This problem is, in fact, a kind of inverse 
of the previous one: now we need to convert the models with characteristic Toda-type 
representation (1.3) into KP-related form (1.2), which is peculiar for GKM. As usually, 
such conversion is provided by Miwa transformation (1.1), moreover the resulting {<fti(m)} 
in (1.2) are also orthogonal polynomials (up to trivial factor). As a particular result, we 
explicitly prove the equivalence of discrete Hermitean one-matrix model and GKM with 
V(X) ~ X 2 and the zero-time n is identified with the size of the matrix in the discrete 
model. 

The fact that discrete matrix models can be naturally embedded into GKM, allows 
one to formulate and study the "double" -scaling continuum limits as internal problem 
of GKM, where it becomes just a question about asymptotic formulas for families of 
integrals. We present some naive results about such interpretation of continuum limits, 
but its explicit relation to conventional Kazakov's procedure (as described in full details 
in [7]) remains still a bit obscure. 

2 Partition function of GKM as a Toda-lattice r- 
function 

5 The particular model (with i_ p = but n ^ 0) has been proposed and discussed to some extent in 
a very recent paper [5], which stimulated us to complete this research. Despite our disagreement with 
the main claim of that paper (that GKM with V(X) ~ X 2 and n > describes c = 1 string model) 
we shall refer in appropriate place to the proper Ward identity (in fact, the system of discrete Virasoro 
constraints [6]), which is derived there. 
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2.1 Interrelation between KP and Toda- lattice r-functions 

The purpose of this subsection is to describe without any special reference to GKM the 
explicit relation between KP-like (in Miwa variables), 

r^[t p ] = A(m) , (2-1) 

and To da- like, 



T~n[t—pitp] — det Hi-\- n j+ n [t_p, t p \, (2-2) 

representations of r-functions, where 



A{m) = Y[{m l -m j ), (2.3) 



( m ) =m i - 1 (l + 0(-)), (2.4) 

7TZ 



* P = -£™r p ,p>o, (2.5) 



c^/ctt p = if^, j > p > 0, (2.6) 



and 



dHij/dt.p = Hi_ pJ , i>p>0. (2.7) 

(the origin of these formulas, concerning with Toda hierarchy, is explained in Appendix 
A). 

Relation between (2.1) and (2.2) is formulated in terms of the Schur polynomials, 
which are defined by: 

V[z\t p ] = exp{]T V^} = 5>*P fc [*], (2.8) 

p>0 

e.g. P_ n = for any n > 0; P [t] = 1; Pi[t] = ti; P 2 [t] = t 2 + \t\\ P 3 [t] = t 3 + t 2 h + ±t? 
etc. The crucial property of Schur polynomials is: 
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dP k /dt p = P k ^ p (2.9) 

(this is just because dV/ dt p = z p V). This feature allows one to express all the dependence 
on time-variables of Hij[t], which satisfies eqs.(2.6) and (2.7), through Schur polynomials 
(see (A. 26) in Appendix)^: 

Hij[t-p,t p ] = Pi^k[t-p]TkiPi+j[tp\, (2-10) 

h-<i 

where T k i = H k i[0,0] is already a t- independent matrix. Note that is defined by 
eq.(2.10) for all (positive or negative) integer values of Somewhat different, the 
Grassmannian point, {<f)i(m)}, in KP-formula (2.1) is a priori defined with i > 0. 

Therefore, we begin our consideration from the case, when all n = i_ p = 0, then look 
what happens if n > 0, discuss the continuation to negative values of n and introduction 
of i_ p - variables. At the end of the subsection we discuss the conditions for forced and/or 
Toda- chain reductions to occur. 

Given the system of basic vectors 4>i(m) for i > 0, we put by definition 

H^t.p = 0,t p ] = <f ^{z)z-'V[z\t p }dz, i > 0. (2.11) 

The integration contour is around zero and it can be deformed to encircle infinity and the 
singularities of V[z], if any. If we just substitute the definition (2.8) of V[z) into (2.11), 
we get (2.10) with Pk-i[t- p = 0] = 5 k i and 

T kl = I (f) k {z)z l dz. (2.12) 

Jz^O 

In order to prove the identity between (2.1) and (2.2) under the condition (2.5), note that 
(2.5) implies that 

r (-) fc A k (m) 
(z - m k ) A(m) ' 

where A k (m) = n i m i ~ m j)- Note that eigenvalues m k = oo do not contribute to 

V[z\tp\. If there are exactly iV finite eigenvalues m k ^ oo, then the point z = oo does not 
6 We manifestly write down the limits of this sum only for convenience as they are given automatically 
by properties of Schur polynomials. 



V[z\t p 



det M 
det(M-Iz) 



n 



(mi - z) 



n 



mi 
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contribute to the integral (2.11), at least, for % < N. It picks up contributions only from 
the poles of V[z\t p ] at the points m k : 



Hij\t—p 0, tp] 



4>i{z)z J V[z\t p ]dz = — — y2^(-) 4>i{m k )- 



The sum at the r.h.s. has a form of a matrix product and we conclude that 



det Hij = det <pi{m k ) ■ J] 



lit ^ 

A(m) 



-) k A k (m) 



det—.. 

mi 



The last determinant at the r.h.s. is equal to A(l/m) = (— ) N( - N 1 ^ 2 A(m) ■ Y\ k m 



Note also that J | 



Afc(m) 
A(m) 



A(m) 2 , and gathering all this together, we see that the 



r.h.s. of (2.12) is indeed equal to 



A(m) 



(2.13) 



as required. 

Proceed now to introducing of zero- and negative-time variables. The zero-time n 
arises just as the simultaneous shifts of indices i and j of : — > H i+n j +n , see (2.2). 
We can use eq.(2.11) to write: 



Hi 



i+n,j+n 



[o,t P ] 



<p\ n \z)z~' J V[z\dz 



(2.14) 



with 



<t>r\z)=z- n <t> i+n {zy. (2.15) 

This exhausts the problem of restoring the n-dependence for positive integer values of 

n. As to n < 0, eq.(2.15) can be used only if original set {<pi(z)} is enlarged to include 

4>i{z) = z l_1 (l + O(-)) with negative i. Such extension is not limited by any additional 

constraints and is not unique: this is exactly all the new information, introduced when 

7 Let us note that the relations = H i+n j +n and cj>^ = z~ n (j) i+n are correct only for the special 
choice of H and <j> (remind that they are defined up to triangle transformation). Moreover, the choice 
of H fixes <p unambigously due to eq.(2.11): H i+n j +n = (<j>^ z~^) = {4$^ z~ k z~i) for any j, i.e. all 
moments of functions <f>^ and (f)^^ 1 z~ k are coincide. This property being evidently correct for GKM 
(f)[ n \ singles out the latter. 
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proceeding from KP to Toda-lattice hierarchy^. If this extension is chosen, Hij\t- P = 0, t p ] 
and the entire matrix Tjy with any integer (positive or negative) i,j,k,l are defined. 
From the point of view of Grassmannian the integer-valued zero-time n labels connected 
components of Grassmannian, consisting of the vector sets {<pi +n (z), i > 0}. 

As for negative-times, as soon as T ki is defined, they are introduced with the help of 
(2.10) and 



k<i 



<f <t>f-*> n \z)z-iV[-\t„ p \P[z\t p ]dz, 



(2.16) 



with 



dt-rri (z) 



' z 



(2.17) 



: " exp \ - t-pZ P \ Yl Pk[t- P }(f) l+ n-k(z) 

p>0 



The role of the exponential prefactor in (2.17) is to guarantee the proper asymptotic 
behaviour 

0f -^ n \ z ) = + 0{-)}. (2.18) 

z 

Note that the quantities defined by eqs.(2.16) and (2.17) depend crucially on 0j with i < 
when "zero-time" n is negative. 

Eqs.(2.16) and (2.17) provide a complete description of the interrelation between KP 
and Toda-lattice hierarchies. Given a KP r-function in the form of (2.1) (i.e. in Miwa 
coordinates (2.5)), it can be interpreted as the Toda-lattice r-function at the vanishing 
values of n and {t- p }. When zero-time n is introduced, it corresponds to discrete shifts 
8 In the fermionic language of Appendix A this means that the form of KP r-function allows one to 
restore the element g of GL(oo), which defines the point of Grassmannian, only up to a part, which 
cancels the vacuum \n) (n = for conventional KP). Therefore, the choice of the negative n vacuum gives 
some additional information about g in compare with n — (the same information can be, certainly, 
obtained from the form of the negative time dependence). Thus, generally speaking, there are plenty 
ways to continue g to "negative part" . But it is not the case if one respects particular reduction. For 
example, the condition of Toda chain reduction unambiguously fixes the continuation to negative n [4] . 
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along disconnected components of Grassmannian, while evolution along negative-times 
t_ p is a smooth movement of a point in Grassmannian, explicitly described by (2.17). (In 
other words, at any given values of n and {t- p } the Toda lattice r-function (2.2) as a 
function of positive-times {t p } can be considered as KP r-function, but the associated 
component and particular point of Grassmannian are different for different n and {t_ p }.) 

Further simplification of (2.17) can be achieved only for some particular choices of the 
basis {(pi(z), i e Z}. A drastic simplification arises for basises, relevant for matrix models 
(i.e. consistent with string equations). Then <j>i(z) are essentially of the form (x l ~ 1 ) z (with 
certain linear averaging operation < . . . > z ) (like contour integral representation)! and 
introducing of n and {t_ p } may be naturally described as a change of the "measure" from 
< ... > z to 



Moreover, this formula allows one to consider n as continuous rather than discrete param- 
eter. This makes possible an "analytic continuation" in n and, thus, implies a "natural" 
definition of fa with negative i. Before we proceed to a more detailed discussion of this 
situation in the next section, let us comment briefly on two important reductions of the 
Toda-lattice hierarchy. 

The first reduction which is of importance as it is just the case in matrix models is 
already mentioned forced hierarchy. It was firstly introduced in [8,4] for the Toda chain 
but can be easily extended to Toda lattice case. The most manifest way to give this 
reduction is to constrain the element g to give some subspace in Grassmannian. This is 
done in Appendix B, but here we would like to say some words in our previous framework. 

That is, instead of half-infinite determinant in (2.2) let us consider finite determinant 
of the size nxn for r n : 



9 This interpretation of string equation seems to be a very promissing view on its meaning. Details 
are, however, beyond the scope of this paper. 




(2.19) 



r n = det Hij. 

ry~i V m J 



nxn 



(2.20) 
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It just gives a r-function of forced hierarchy. Moreover, due to identities (2.6)-(2.7) it 
can be rewritten asQ 

r n = det d^dP^H (2.21) 

nxn 

with d = d/dti, d = d/dt-i, H = H\ t \ and H is constrained to satisfy dH/dt p = d p H, 
dH/ dt^p = d p H. This r-function can be produced in formalism of orthogonal polynomials 
implying polynomial Baker- Akhiezer function as we shall see in the sect. 3.1. It explains 
why matrix models correspond to forced hierarchies. 

For all these it is crucially to work in the sector with positive n. The problem of con- 
tinuation of forced hierarchies to negative zero-time we have already discussed previously 
[4]- 

Another important reduction from Toda lattice is Toda chain (see, for example, [10]). 
It can be easily written both in terms of element g ([g, Jk + J-k] = - see (A. 29) and 
comments there) and in determinant form. Latter one merely implies the symmetry 
property: 

[H , A + A" 1 ] = 0, (2.22) 

where A is shift matrix A^- = This condition leads to r-function of Toda chain 

hierarchy (proper rescaled by exponential of bilinear form of times) which depends only 
on the sum of positive and negative times t p + but not on their difference (one 

can consider this as defining property of Toda chain hierarchy). Let us remark that 

one possible solution to constraint (2.22) is matrix Hij = 7ii-J^. We can combine both 

10 Let us point out that this kind of solutions to Toda lattice hierarchy was firstly invented by Leznov 

and Saveliev [9]. Let us also note that solutions of Wronskian type [10] rather like forced ones turn out 

to be absolutely different in their properties. 

11 Let us emphasize that usual notations correspond to the reduction to iait-function independent of 

the sum of times. But due to non-standard sign of the first exponential in (1.7) which simplifies a lot of 

formulas in this paper, the reduction (2.25) corresponds to r-function which depends only on the sum of 

times. 

12 It should not be certainly independent of the difference of times t p — t—p, but one can through out 
negative times as the final answer for complete objects like r-function should be really independent of 
this difference. 
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reductions to reproduce forced Toda chain hierarchy. In this case one can easily transform 
Hi-j to matrix H-i+j by permutations of columns what does not effect to the determinant. 
This matrix just corresponds to one-matrix model case [4,6] (see also sect. 3.1). Thus, we 
consider again the determinant of size nxn, which now can be represent in the form: 

r n = det d l+j H (2.23) 

nxn 

where d = d/dti, dH/dt p = d p H and we canceled negative times (see footnote 12). Like 
the Toda lattice case, the forced Toda chain is unambiguously continuable to negative 
values of zero-time leading to (1.5)0. 

2.2 GKM in the context of Toda lattice hierarchy 

The purpose of this subsection is to introduce zero- and negative-time variables into GKM 
in such a way that its partition function becomes a r-function of the Toda lattice hierarchy. 
Let us remind the definition of GKM at n = t_ p = 0. It was introduced in [1] by the 
following matrix integral: 

r e -Tr{V(X+M)-V(M)-V'{M)X} jj£ 

z ^ m = - — j^Fw^dx — (2 - 24) 

with NxN Hermitean matrices X, M and V 2 {X, M) = lim e ^ e" 2 [V(M + eX) - V(M) - 
eV'(M)X]. As explained in details in [1], the r.h.s. of (2.24) may be represented in the 
form of (2.1) and thus Z{v}[M] is the KP r-function, which has no explicit dependence 
on N (N is just the number of finite eigenvalues of the matrix M, when others can 
be considered as infinitely large). The relevant set of functions {<fti(m)} - the point in 
Grassmannian - is given by the following integral formula: 

(2.25) 

ee 8 {m) J dx x i-i e -v(*)+xV(m) _ 

13 This continuation is unambigous only with taking into account the negative times. With cancelled 
negative times, i.e. in KP case, one can obtain plenty different continuations like, for example, CKP 
r-function which satisfies T n [tk] = r_ n [(— ) k tk] [4]. 
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(the integral in (2.25) is a contour integral, i.e. x not a matrix). Dependence of n and t- p 
is now introduced by the rule (2.19)0: 



,P>0 



x p — m p 




fdx x n+i-l e -V(x)+xV'(ra) exp J ^ £ 



py x y — m 



where 



X p 

V{X)=V{X)-n\ogX-Y,t- P X~ P = E s v— > ( 2 - 27 ) 

p>0 p=— oo P 

with 

t p = ^[TrM-P+(p-l) Sp ], 

t =n = -s , ( 2 - 28 ) 

As to original potential V(m) it can be identified with V + . We repeat that from the point 
of view of GKM to = —so does not need to be integer (though it is desirable to avoid a too 
complicated analytical structure of Z[M\ - see (2.29) below). From (2.26) we immediately 
conclude (just repeating the arguments of [1] in the opposite direction) that the partition 
function of GKM, involving zero- and negative-times (and automatically being a Toda 
lattice r-function), is justQ 

14 Let us point out that the exponential of negative powers in normalization does not essentially ef- 
fect to the KP r-function as it reduces to trivial exponential of bilinear form of times in front of r- 
function and corresponds to the freedom in its definition [1]. Indeed, t ~ det {exp[^ fc a^zl \4>i(zjj} ~ 
EL ex PE a kzf k ] det (f>i(zj) ~ exp[^ fcafeffe] det ^)j(mj). This factor, certainly, effect to the string equa- 
tion, however, results into trivial bilinear form of times (see subsection 2.3, eq.(2.37)). 

15 All the integrals are certainly defined by analytic continuation and are not quite unambigous because 

of Stokes phenomenon etc. 
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frIY „-TrV(X)+TrV!(M)X 
Z fin [M] = e TrV(M)-TrMV>(M)l^_^ _ 16 (2 _ 2g) 

Occurrence of GKM in the somewhat unexpected form of eq.(2.29) (involving projection 
V — > V + ) makes the whole subject even more intriguing than it was before. Our next 
purpose is to describe the associated modification of the second crucial ingredient of GKM: 
the string equation (the first ingredient is integrability). 

2.3 Generalized string equation 

Two ways to derive string equation and W- (in particular, Virasoro) constraints were 

suggested in [1]. One immediately leads to the string equation Z = and, as a 

consequence of it together with an additional information about Z (that it is a r-function 

of a specifically reduced KP hierarchy), one reproduce the set of all other ^-constraints 

appropriate for the given reduction [11]. Another way [12] starts with exhaustive matrix 

Ward identity [13]Q which is transformed then into the set of appropriate ly-constraints, 

without addressing to properties of Z at all. Below we concentrate on the first approach 

which is much simpler in general situation. 

A r-function is parameterized by a point of Grassmannian. The role of string equation 

is to specify this dependence for the case of r-function given by the GKM partition 

function. Unlike the Toda case, if we look at Z^yy as at KP r-function, the point of 

Grassmannian depends on entire potential V(X), i.e. not only on the polynomial piece 

V(X) = V + (X) but as well as on its negative and zero times (being included in definition 

of the point of Grassmannian in KP framework). 

16 Note also that Z^ v y [M] ^ Z{v} [M] because of occurrence of V+ at some places. One more model can 
be denned by the rule .F^j[V|.(M)] = F^y[V' (M )] (where JF's denote only integrals in the numerators 
of (2.24) and (2.29)). Then M is defined by the relation V{{M){= V'(M)) = V'(M), and Z {V} [M] is a 

KP r-function of somewhat different time variables t p = — \TrM~ p + (p — !)§„], where s p are defined to 

V 

be the coefficients of V expansion in powers of M. In the simplest case of V(X) = -X 2 -nlogX the 
definition of M is M — M — nM _1 . This is a rather familiar relation in Toda chain theory connecting 
its spectral parameter with that of KP hierarchy. We emphasize that while both Zj-yj[A/] and Z^y^[M] 
are KP r-functions, only Z^yy is also a Toda lattice r-function (i.e. possesses a proper dependence of to 
and i-p). 
17 See also [14] and [15]. 
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In this subsection we shall shortly repeat the derivation of string equation^ addressing 
reader for more details to [1]. We shall demonstrate that this derivation crucially uses 
only the following information: 

1) given asymptotics of <j>i(ji) ~ A** -1 ; 

2) the manifest form of normalization factor depending on M in GKM integral; 

3) the manifest form of linear term in X in exponential in integrand (2.25) (it is 
necessary in order to get correct expression (2.1) for Z implying Z to be r-function). 

The main idea is to consider the following derivative of the GKM partition function 

7 _ detUmj) _ det s( mj )U™j) ^ (2 24) _ (2 25)): 



A(m) Aim) 



and rewrite it as 



y Tr 1 dt p dlogZ _ ^ 1 1 dlog Z 

^ V"(M) 8M dt p ^ V"(M) Mv+i dt p K - J 

On the other hand, the derivative (2.30) is equal to two pieces: the first one originates 

from the derivative of factors s(m) and A(m) _1 (here we use the information of point 2) 

and is equal to: 



2 V"( mi )V"( mj ) mi 



m 



j 



The remaining second piece can be transformed to derivative over t\ essentially using the 
correct asymptotics of (piijn) (point 1): 



Tr i 77^777V7777 lo § det H m i) \ = Z ■ ( 2 - 33 ) 



V»(M)dM l ° g det ^K')} = ^ 
Now let us modify GKM partition function by introducing of negative- and zero-time 
variables, in accordance with (2.29). Then the only change (one should certainly to 

replace all V by V + ) of string equation originates from the additional piece in (2.32): 

18 In this section we really deal with L_i-constraint, the derivative of which over the first time variable 
is usually named string equation. Nevertheless, for the brevity, we call this Virasoro constarint as string 
equation too. 
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Tr^i — dV ~?f^ log Z = -t Tr-„ — + y>_ p 7V^ . (2.34) 

Thus, we finally obtain the string equation in the form: 



L_!T = 0, 

L {v+} - sr- T {v + } + ^_ _ T {v+} ,sr- T {v + } f , 



with 



71 { ^ } = Tr^ . (2.36) 

p V^(M)Mp 

Let us note again (see footnote 14) that it is not necessary to include negative times 
into normalization factor s(m). In this case, the string equation is not at all modified by 
negative times. Moreover, we can multiply GKM integrand by exponential of arbitrary 
series J2-oo a pX p with the only restriction k < (the leading degree of V(X)) to preserve 
the property 1 of correct asymptotics of 4>i- One can easily see from our previous dis- 
cussion that this does not disturb three essential properties above and the form of the 
string equation is conserved. Thus, in KP framework we have plenty deformations of 
starting point of Grassmannian not changing the string equation. These modifications 
does certainly change the reduction condition, therefore, it emphasize the importance of 
fixing the reduction to define uniquely the r-function of (reduced) KP hierarchy by the 
string equation, in spirit of [11]. 

On the other hand, one can merely consider this as an extension of standard GKM 
viewpoint. Nevertheless, usual GKM representation is singled out by the possibility to do 
smooth transition to potentials of the other degrees corresponding to the double scaling 
continuum limit of the other multi-matrix models. In this sense, one can look at the string 
equation as giving some sort of universality class. 

Now let us return to the string equation (2.35). For particular potential V + (X) = 
X K+1 /(K + 1) , T p+ i — (p + K)t p+K and L^-constraint has a form: 
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L { -V = Y,(P + K )t P+ K^- + K^-+ 1 - £ at a bt b (2.37) 

p>0 Ul P ° l1 Z a+b=K 

with a, b being any (positive and negative) integers and at a should be substituted by n 

for a = 0. For all t- p = the term d/dti can be interpreted as resulting from the shift 
K 

t P ^t p + — - - <)p,A-+i = V K also K — 1, 

L { -1 = E(P + !)W J" + w *i ( 2 - 38 ) 

p>0 OT P 

and this is similar to the string equation of discrete Hermitean one-matrix model, n being 
the size of the matrix. We shall return to this point a bit later. Here only note that, 
while in [5,6,16] it was proposed to interpret n as d/dt , we see now it is more natural 
identify n with t itself, in the sense that it is n that plays the role of zero-time in Toda 
hierarchies (this was mentioned but not emphasized in [6]). 

Somewhat alternative approach to the derivation of string equation and associated 
tower of WMike constraints could be to begin with the Ward identity for the integral 

F(A) = JdX e -TrV(X)+TrAX ? ^gg) 

which results from the invariance under a shift X — > X + e (e being a matrix) of the 
integration variable: 

(Tr e{V'{X) - A}) = 0. (2.40) 

Usually we can represent the positive powers of X coming from V+(X) in the integrand 
(before averaging) in (2.40) by action of derivative operator V+(d/dA t r) to already aver- 
aged expression. This is naively impossible for negative powers of X in V!_(X). However, 
be there only finite number of negative powers (i.e. if V!_(X) is a polynomial in X -1 ), one 
could take additional A-derivatives of equation (2.40) in order to obtain more factors of 
X and eliminate its negative powers. This trick was applied in [5] in the particular case 
of all t_ p = with only non-zero n. In this way they obtained the identity of the form: 



d 



tr 



d 

+ 



A + nlF(A) = (2.41) 



dA tr J 
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which, if rewritten in terms of Z and in the limit of N — > oo (see details im [12] and [5]) 
turns into a set of Virasoro constraints 



L P Z = 0, p > -1 

with 



d ™ _1 d 2 d 
L p = H kt kjf f + E + 2n— for p > 1, 



d 

L = E o 1" ^ 2 an d -^-i as m (2-38), 

fc=i 



(2.42) 



which are identified with Virasoro constraints for Hermitean one-matrix model. If some 
finite number of negative times are non-vanishing, one get instead a set of VF-constraints 
[17]. 



3 Discrete models in the form of GKM 

Since we devote this paper to discussion of Toda lattice hierarchies in the context of 
matrix models, we can not avoid touching the main conclusion of ref.[6] that all the 
discrete matrix models do correspond to particular cases of Toda hierarchies. In the 
simplest case of Hermitean one-matrix model one gets a Toda chain, other multi-matrix 
models correspond to other reductions of the Toda lattice hierarchy. Moreover, all discrete 
matrix models fall into the class of forced hierarchies [4]. 

3.1 Discrete models as forced Toda-lattice r- functions 

In this subsection we show how partition functions of various matrix models can be 
rewritten in the KP-like form of eq.(2.1). 

The first step in this direction is to reproduce them in the Toda like form of (2.2). This 
statement was already discussed in [6] and [4], but we present it more explicitly below. 
The second step is to notice that ify which arises from discrete models are representable 
as averages, namely, finite dimensional integrals, or "matrices of moments" : 



H tJ [t^ p ,t p ] = (xyV[x\t p ]V[y\t^ p ]). 
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The third step is to perform a Miwa transformation of times t p with t- p fixed, so that H^s 
become averages of polynomial functions of X. Then det may be transformed with the 
help of orthogonal polynomials technique. These polynomials are, however, orthogonal 
with respect to the simple measure, defined by potential V(H) = V(H) —HV'(H), which 
depends only on the s p - variables and not on the times t p . 

The main result of all these calculations is that partition functions arise in the form 
(2.1), with 4>i{z) being proportional to orthogonal polynomials. 

First, we shall briefly repeat what is known about discrete matrix models. In the case 
of the Hermitean one-matrix model the n x n matrix integral 

ZW = {Vohwnl}- 1 J DH exp{-]T^p/P} = 

= (niy 1 fUdhAWeM-Y,^} (3 ' 1} 

i i,k 

is taken by orthogonal polynomials for the arbitrary potential W(h) = J2tkh k 

< pf\ P V >= / P?\h)pf\h)e- W ^dh = tf ye wW (3.2) 

and equals 

Z (D = jj eW W =T (i) (0 (3.3) 

which is the r-function of the forced Toda chain hierarchy (being defined for negative 
n by (1.5)). In eq.(3.1) Volu^ denotes the volume of unitary group U(n) (see (3.41) 
for explicit expression). On the other hand, it follows from (3.2) and the definition of 
orthogonal polynomials 

^ (1) W = E«S ) ^' 

j ^ (3.4) 
a? = 1 

that 

diag(e^) = A^C^A^ T , (3.5) 
where = ||a^||, A T - transponed matrix, and is so called matrix of moments 
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Thus [4], 



C \f = J h^e~ w ^dh. 



(3.6) 



T W(t) = det[diag(e^)] = det A^C W A^ T = det C {1) . (3.7) 

More or less the same is true for all discrete matrix models, only the form of A (which is 
unessential in (3.7)) and scalar product in the definition of C depend on the model. In the 
case of two Hermitean matrices and two potentials W(h) = J2tkh k and W(h) = J2ikh k 
we have 

= {Volu^nl}- 1 J DH DHe-Splww+ww-Hfi] = 
= N)~ 1 |nrfM^A(/i)A(/i)exp{-^[Vy(/ il ) +W(Jn) -hihi}} = (3 g) 

n-l 
i=0 

where 

e ^m 5 .. =< P (2) ) pW >= J P p(h)pP(h)e hh - w ^- W ^dhdh (3.9) 

and 

ri 2 \t,i) = det[diag(e^)] = det A^C^A^ T = det C (2) (3.10) 
for matrix of moments given by 

C \f = J VWe h ~ h - w ^- W ^dhdh . (3.11) 

In generic i-T-matrix situation we can also use the formulas (3.10) and (3.11) if we take 
care of the integration measure. Namely, 
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K-l K 



Zf» = {Vol u{n) n\}- 1 j DH DH JJ DH® exp{-Sp{W(H) + W{H) + J2V (l) {H (l) )- 

1=2 2 

- J2 ciH^HV+V - HH® - H^H}} = 
= (n!)" 1 / n rf M^ n 1 ^f ) A(/i)A(/i)exp{-^[^(/ il ) + ^(^) + E^^ )- 

J i 1=2 i 2 

n-l 

= jj eVi (t,tiW( l ) lCl ) =r {Ar iW (0 )Cl}[fjf -| ^ 

i=0 

(3.12) 

where we distinguished if^ = H and = H among other integration variables, 

because the role of Toda time variables is played by = t and = i, while other times 
t®; I = 2, . . . , K — 1 (or potentials W®) and the set of {q} are considered as parameters. 
(The evolution along these frozen times can be described by the multi- component Toda 
lattice hierarchy.) 

The orthogonal polynomials now look like 

e <Pi(t,t;{W},{c}) 5 .. =< p.W ; pW > = 

(3-13) 

= Jp t {K '(k^W^lM] , 

where 

d^ K \h,h) = dhdhe- s ^ w ^ +W ^ [Y[ dh® exp{-Sp{xV«(/i«)- 

7 2 2 (3.14) 



Again 



T {K; { w},{c}} [t ^ = det[diag(e^(*'^ H/ >'W))] = 

(3.15) 

= det aWcWAWt = det C W 



and the matrix of moments is 



= Jh l Wd^ w ^ c »[h,h} 
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(3.16) 



In the limit K — > oo this should correspond to the infinite-dimensional matrix model, 
whose first critical point has to describe (at least, formally) c = 1 conformal matter cou- 
pled to gravity [18]. Instead of the infinite sum we introduce the "quantum mechanical" 
integral and (3.12), (3.14) can be rewritten as 

w ~ f°^nmQert-sp£dt{w[H®] + (f ) 2 } } = 

= / Ud^\hi,hi)A{h)A(h) = (3.17) 

n-1 

= JJ e w[MiW(0] =T W)}[ M ] 
i=0 



with 



d^\h,h) = dhdhj h ^ =h HDh(t)exp{-Spj o ddw[h(0} + i^J |}. (3.18) 

So far we considered the partition functions of various discrete matrix models as functions 
of Toda time- variables. Now we shall demonstrate that if one passes from times to Miwa 
variables 



1 1 N t> — 1 

t„ = -TrM~ p = -Y m- p + Z s p (3.19) 

P V~{ P 



(j) > 0) and/or 



1 1 N v-1 
t = i = -TrA~ p = - V \~ p + P - s p (3.20) 

P Pf^i P 

(note that N - the size of matrices M and A has nothing to do with n - the size of 

matrices H, being integrated over in (3.1), (3.8) and (3.12)) then the partition functions 

of discrete models {t^} acquire the form of KP r-function (2.1). Indeed, 
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i i,k 



hi 



= (n\)- 1 / HdhA 2 (h)e- v ^H(l -^-) = 



71! 



i-l 



A(h,m) 
A(m) " 



x det P^Uhj) det 

raxn (n+AT) x (ra+JV) 



(0 



k — 1 



(3.21) 



where i, j = 1, . . . , n; a, b — 1, . . . , N; V(h) = ^ ' s kh k and {P^\h)} are correspond- 
ing orthogonal polynomials with respect to deformed measure e~ v dh : 



< pW,pW >= f pW(h)pV(h)e-*Wdh = V ftW 



so that 



(3.22) 



Pl l \h) = h l + 0{h i - 1 ). 
Computing determinants in (3.21) and using orthogonality condition (3.22) one obtains 



r^Hs] = nVA-^m) det P^ a _iK)Il e 



(i) 



3 



n 



NxN 



det (a6) ^ 1,n) (m 6 ) _ (1) det (ab) ^ 1>n >(m 6 ) 

_ — T n P°I S J X _ 



(3.23) 



A(m) n A(m) 

i.e. the r-function of the discrete Hermitean one-matrix model acquires the form of eq.(2.1) 
with 



m ) = m n P^+ ) _ 1 (m) 



(3.24) 



Below we shall see that (3.24) is natural representation for all discrete matrix models. 

Let us remark that the expressions (3.21)-(3.24) does not depend on the quantity N. 
It means that one can consider eq.(3.21) with N = 1 and reproduce well-known integral 
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representation for orthogonal polynomials [19] (the serious drawback of this representation 
is that its manifest form, and, moreover, the number of integrations depend on the degree 
of polynomial): 



In the case of two matrices instead of (3.21) one gets 



rf (t,t) = (n!)- 1 / \{dhdhA{h)A{h)e W {-Y.\ W ^) + W ^ l ) - hJa]} = 
J i i 

= (nl)- 1 [ 1] dh t dh t A(h)A(h) exp{- ]T[F(/0 + Hh) - M]} W 1 ~ ~ ) = 

J i i i,a m « 



A(h,m) 
A(m) 

(nVt^Hm^A-^m) J J[ dhidhi exp{- ^2\V(hi) + Vfc) - M]}x 

i 



-(2) - 

xdetP.-iiL det 

nxn ■>' (n+N)x(n+N) 



n 



3 V»(s,*) 



det (a6) 



AM 



where 



(3.25) 



0i 2 ^(m)=m-"Pi 2 + U(m), 



(3.26) 



< p/ 2 ), jf >= / pP(h)pf\h)e h ~ h - V ^- V ^dhdh = e^ s '% . (3.27) 

Note that, using symmetry of the two-matrix model, one can obtain the formula (3.25) 
using Miwa-transformed form of the other set of times (3.20), though its proper interpre- 
tation is not yet completely clear. In this case: 



n 



det (flt) ffi">(A 6 ) 
A(A) 



(3.28) 
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with 



(3.29) 



where 



< Q?\Q? >= J Ql 2 \h)Qf\h)e h ~ h -^-^dhdh = e^5„ 



(3.30) 



for V{h) = Y J 1 ^s k h k . 

Let us point out that formulas (3.25), (3.28) look especially nice in the particular case 

of asymmetric two-matrix model [17,20] with one of the potential being a finite polynomial 

M k - 1 ~ 
of fixed degree, say W(h) = ^ — ; — s k h k . Then it is natural to take V(h) = W(h) and 

k=i k 

(3.25) turns to be 



where 



n 



3 Vi(s>s) 



det {ab) ^ n ' w \m b ) 
A(m) 



(3.31) 



(3.32) 



<P^ W \P^ W) >-- 



(3.33) 



= J P^ W \h)P^ W \h)e hh - w ^- w ^dhdh = 
As we have seen above orthogonal polynomials technique without serious changes can be 
applied to generic if-matrix model. For the analogues of (3.25), (3.28) it gives 



n 



L i 



det (a6) <Pi K ' n >(m b 



A(m) 

det (a6) ^ n \\ b ) 
A(A) 



(3.34) 



where 



<P { a K > n \m)=m- n P£Um) , 



(K) 



(3.35) 



26 



< Pf \ ~P { P >= I P[ K \h)P^(h)e h ~ h -^- V ^d^ K \h : h) = , (3.36) 

&*.»>(A) = A-»CL(A), (3-37) 

< Q?\ ~Q { P >= J Qf\h)Qf\^ hh - v W-^ h W K) (h, h) = . (3.38) 
3.2 Hermitean one-matrix model 

In the previous section we demonstrated the manifest connection between r-functions 
of discrete matrix models and KP r-function in Miwa variables. However, in order to 
be representable as GKM they still need to arise in a somewhat specific form. Namely, 
components of the vector {0;(m)} should possess a representation as "averages", 

Um) = (x 1 - 1 ) . 

\ I m 

Since in the study of discrete matrix models (f>i(m) arise as orthogonal polynomials 
P n+i (m), what is necessary is a kind of integral representation of these polynomials, 
with i-dependence coming only from the x' -1 -factor in the integrand. It is an interesting 
problem to find out such kind of representation for various discrete models, but it is easily 
available only whenever orthogonal polynomials are associated with the Gaussian mea- 
sure: the relevant Hermit polynomials are known to possess integral representation, which 
is exactly of the form which we need. We saw that the choice of the measure is rather arbi- 
trary in the study of discrete matrix models, as soon as the complicated time-dependence 
is eliminated from the measure with the help of Miwa transformation. Therefore, it may 
be possible to make this measure Gaussian. The obstacle can arise either in multi-matrix 
case, if we do not want to apply Miwa transformation to "intermediate" times, or in the 
interesting version of two-matrix model with one potential fixed [17,20], or in the models, 
involving non-Hermitean matrices. Each of these situations deserves special analysis. No 
more details are necessary in the simplest case of one-matrix Hermitean model, to be 
discussed below for illustrative purposes. 
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The statement to be discussed is that Hermitean one-matrix model with the matrix 
size n is equivalent to GKM with V(X) = X 2 /2 — n logX. We can come to this conclusion 
in various ways. 

First of all, we can look at the complete set of Virasoro constraints for this model, 
following from the Ward-identity (2.41). It has been derived in [5] and shown to coincide 
with the set of Virasoro constraints for Hermitean one-matrix model, as discovered in 
[6,16,21]. 

Second, we can prove that such GKM corresponds to a forced Toda- chain hierarchy. 
This statement, if combined with the string equation (2.38), would also lead to the same 
conclusion. 

Third, we can just verify the explicit identity between Gaussian matrix integrals, 



JdH nxn det(I-H/M)e- s P H2 / 2 J dX NxN det(/ - iXjMYe^l 2 

fdH nxn eSpHV2 f dX NxN e~TrXV2 ■ ^) 

Note that the size of matrix in the l.h.s. is n x n and in the r.h.s. is N x N, and these 
parameters are absolutely independent. This identity is indeed true for any n and N, as 
follows from the reasoning of the previous subsection 3.1, and integral representation of 
Hermit polynomials (fij(^) = Hej(ifi). 

Remarkably, we could use this explicit proof (of eq. (3.39)) as an manifest verification of 
the commonly accepted belief that either («) the string equation (L_i-constraint) plus the 
fact that it is imposed on appropriately reduced Toda lattice r-function or (ii) the entire 
bunch of Virasoro (or W-, if necessary) constraints (without any a priori information 
about integrable structure) defines the partition function uniquely (and, in particular, 
predetermines the validness of all other H^-constraints, suitable to the given reduction, 
in the case (i), or guarantees that partition function is, in fact, the appropriate Toda 
lattice r-function in the case (ii)). This belief is, of course, implicit in the first two of 
the above "derivations" of identity between Hermitean one-matrix model and GKM with 
V(X) = X 2 /2-n\ogX. 

In the remaining part of this subsection we shall prove that partition function of such 
GKM is indeed a r-function of forced Toda chain reduction of Toda lattice hierarchy. 

Let us begin with the word "forced" . Conceptually, the relevant hierarchy should be 
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forced just because we deal with discrete matrix models, orthogonal polynomials, and, 
thus, all 4>i(m) for i > should be polynomials. It is, however, somewhat more complicated 
to prove the property (1.5). It even looks a bit intriguing, since partition function of GKM 
(the integral at the r.h.s. of (3.39)) does not seem to vanish for negative n. The resolution 
of the puzzle comes from the study of n-dependence of the coefficient of proportionality. 

To do this, let us accurately restore all normalizations in r-function and use formula 
(3.39): 

_ _ 1 f H „ p -TrV {H) _ (27T)" 2 / 2 J dX NxN det(I - X , ' M)e~ T ^ V 2 

Tn n\Vol u{n) J dUnxne 'n\Vol u(n) ~~ fd x NxN e-W " ' [ " ^ 

where Voljj( n ) denotes the volume of unitary group (it appears due to necessity of inte- 
gration over angular variables), 

n 

Vol~\ n) = (27r)-«(- 1 )/ 2 TT kl . (3.41) 

k=l 

This formula can be obtained, for example, by immediate calculation of matrix Gaussian 
integral I = J dH e~ TrH2 = (27r) n2//2 and comparing it with one calculated in orthog- 
onal polynomials technique: / = Volu( n ) / II dmie~ m ^ A 2 (m) = Volu( n )n\ I\k=i( norm °f 
Hek) = Volu( n )(2ir) n / 2 rife=i k\. Now we can continue it to negative n using the formula 
(see for example [22]): 

if f(p) = f[<t>(i) , then f(-p) = jf 0(— z)" 1 . (3.42) 

i=l i=0 

It gives in our case 

\n\ 

n\Vol u(n) - (27r)«(«- 1 )/ 2 J] r(-<) , (3.43) 

n<0 

what leads to the statement (1.5) due to singularities of r-function at negative integers. 

Our next point is to prove that the GKM integral under consideration is indeed a 
Toda chain r-function, or, what is the same, = Tii-j. We neglect, for a moment, the 
negative-time variables. To calculate, we use the following Ward identity for the average 
with potential V(X) = X 2 /2 (the logarithmic term is absorbed into the subscript of (f>) : 

(x') z = (x'- 1 z) z + (t-l)(x l - 2 ) z . (3.44) 
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As H, 



' <Pi(z) 



and 4>i(z) = (x l ) z , one can write down: 



H 2 j - H ltj _ 1 , 



H 4j - H 3 j-i + 2H 2 j , 
H 5j = H±j_i + 3H 3 j , 

Then, using the admissible triangle transformation 

Hij = Hij , 



i.e. 



H 2j = H 2j , 



H 3j = H 3j - Hij , 



H 4 j = H A j — 2H 2 j , 



H 5 j = H 5 j — 3H 3 j , . . . , 



k 



(3.45) 



with lower-triangle matrix 



-1 1 

-2 1 

-3 1 



V 



/ 



one obtains the property = Hi-j. It is evident that switching on the negative times as 
well as logarithmic term modifies Ward identity (3.44) by contributions of smaller powers 
in x and can be also absorbed by proper lower-triangle transformation. It completes the 
proof. 
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4 Double-scaling continuum limit in terms of GKM 



We demonstrated in the previous sections that the discrete Hermitean one-matrix model is 
equivalent to GKM with V(X) = X 2 /2 + nlogX, while from [1] (see also very instructive 
review [23]) and [12,14,15,24] we know that its double-scaling continuum limit is described 
by GKM with V(X) = X 3 /3. Thus, we conclude that 



lim Z { y } 

d.s. n^oo 1 ' 



7 2 

z {v}- 



(4-1) 



This relation should certainly be understandable just in terms of GKM itself. Moreover, 
since GKM is not sensitive to the size of the matrices N, it needs to become just a simple 
relation between finite dimensional integrals. 

In this section we present some ideas about this relation, without going into details 
concerning its compatibility with conventional description of double-scaling limit [7]. 

To begin with, let us recall that double-scaling continuum limit for the model of 
interest implies that only even times t 2 k = —rTr should remain non-zero, while all 
odd times t 2 k+\ = 0. This obviously implies that the matrix M should be of block form: 



M 



(4.2) 



M 
-M 

and, therefore, the matrix integration variable is also naturally decomposed into block 
form: 



X 




(4.3) 



Then 



Zj {V=X 2 /2-n\ogX} ~ 

1 ( 4 - 4 ) 

= J dxdyd 2 z &et(xy - z± 2 y) n e- Tr ^ 2+x2 / 2+y2 / 2 - MX+My l 

To take the limit n — > oo, one should assume certain scaling behaviour of X, y and Z. 
Moreover, the notion of <iou6/e-scaling limit implies a specific fine tuning of this scaling 
behaviour. So we shall take 
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X = a(i/3I + x), 
y = a(-ipi + y), 

(4.5) 

Z = a(, 
M = aT 1 ^/ + m) 

with some large real a, f3 and 7. If expressed through these variables, the action becomes: 



Tr{\z\ 2 + x 2 /2 + y 2 /2 - mx + My - n \og{xy - z^-zy)} = 

1 2 1 2 

^-Tr{(i(3I + x) 2 + ^-Tr(if3I - yx) 2 + -f 2 \z\ 2 } - Tr(ial + m)(2i/3I + x - y)- 
-nTr log /3V{1 - + f 2 ~ + o(l/ (3))} = 

(4.6) 

= [2ap - /3 2 7 2 - 2n log /^Tr I - 2i(5 Tr m+ (A) 

+i{(3 1 2 -a + -)(Trx-Try) + -( 7 2 - — ){Tr x 2 + Tr y 2 )+ (B) 
P 2 fi 2 

+ (7 2 -^)Tr|C| 2 - (C) 



-Trmx + Tr my + —Tr(x 3 -y 3 ) + (D) 

+0(n/(3*) + 0(\(\ 2 p. (E) 

We want to adjust the scaling behaviour of a, f3 and 7 in such a way that only the terms 
in the line (D) survive. This goal is achieved in several steps. 

The line (A) describes normalization of functional integral, it does not contain x and 
y. Thus, it is not of interest for us at the moment. 

Two terms in the line (B) are eliminated by adjustment of a and 7: 

2 n 2n , A 

^• Q = 7P < 4 ' 7 > 

As we shall see soon, 7 2 = n/(3 2 is large in the limit of n — > 00 . Thus, the term (C) 
implies that the fluctuations of £-field are severely suppressed, and this is what makes 
the terms of the second type in the line (E) negligible. More general, this is the reason 
for the integral to split into a product of two independent integrals leading to the 
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square of partition function in the limit n — > oo (this splitting is evident as, if Z can be 



neglected, the only mixing term logdet 



( X Z 



turns into log Xy = \ogX + logy). 



z y 

Thus, we remain with a single free parameter (5 which can be adjusted so that 



V 



(3 3 

> const as n — > oo {i.e. (3 ~ y/n), (4.8) 

making the terms in the last line (E) vanishing and the third term in the line (D) finite. 

This proves the statement (4.1) in a rather straightforward manner, without addressing 
directly to the complicated matter like Kazakov's change of time variables, reformulation 
of Virasoro constraints and so on [7]. We do not go into more details here, but point 
out one important detail. That is, the possibility to eliminate all the terms of original 
potential V + (X) of degree K = 2 by the contributions coming from expansion of logarithm 
in such a way that the (K+ l)-th power of expansion survives is due to careful fine tuning 
of parameters of original V + (X). This is just the idea involved in the notion of double- 
scaling limit and for higher-degree potentials it should be replaced by ".fT-scaling" limit 
which turns Z{ V=X K_ n \ ogX } into Z 2 v=xK+1 y as n — > oo. 



5 CONCLUSION 

To conclude, let us first remind the main idea of the paper. The proposed Generalized 
Kontsevich Model seems now to be close to a unified theory of all (discrete and continuous) 
matrix models. By introducing of "negative-" and "zero-" time variables it results to be a 
r-function of the Toda lattice hierarchy and allows one to unify all matrix models in this 
framework. In particular, the original Hermitean one-matrix model is, in these terms, the 
GKM with "trivial" potential X 2 (non-trivial models start from original Kontsevich's one 
with X 3 potential) but with non-vanishing zero-time. This must give an opportunity to 
study all subtle questions in this unified language (the continuum limit of discrete models 
among them - the sketch of such treatment was proposed in the sect. 4). Besides, the 
introducing of zero-time variable leads to the notion of "natural basis" in Grassmannian 
with property 4>i n \z) = z~ n (f)f^ n (z). This basis differs from the canonical one (which is 
singled out in the Segal-Wilson construction [27] and corresponds to standard represen- 
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tation of fa through the fermionic averages in spirit of GKM [1]) and coincides with one 
inherited from GKM integral (2.25). 

One of the most important items is the matter of the sect. 2. 3 concerning with (gener- 
alized) string equation. We see now that the string equation does not determine the point 
of infinite-dimensional Grassmannian uniquely - it allows one to deform it in accordance 
with negative-time Toda flows. In this sense, the introducing of the negative times would 
rather correspond to handle-gluing operators being added to action of string model, while 
the positive times to the motion in the space of string models. What is really spoilt by 

these negative times is the particular reduction condition, which is important in order to 

(q - q') 2 

determine a particular string model with c = 1 — 6 . This condition is also spoilt 

qq' 

by perturbation of the pure monomial potentials X K+1 by lower order irrelevant operators 
(in the sense of the universality class defined by the string equation), however the physical 
properties of these "deformed" models should be the same as those of the original one. In 
particular, it means that the choice of normalization of the basis in Grassmannian given 
by GKM <]>&) ~ J dx e -v(x)+v'(z)x could be more flexible. 

We have only touched in this paper the most serious problem of the continuum limit, 
which has been done correctly up to now only for the Hermitean one-matrix model [7] . As 
it has been shown in the sect. 4, this continuum limit in principle can be done by purely 
"field-theory" technique in the framework of GKM, though some questions are still not 
completely clear. It allows us to hope that the proposed in [1] and above formalism 
will shed light on the serious problem of the continuum limit in multi-matrix models, 
among which the asymmetric two-matrix model proposed in [17,20] seems to be the most 
perspective one. 

Finally, we have to repeat that we do not agree with the interpretation of [5] that the 
potential X 2 — n log X describe consistently c = 1 string theory, or, equivalently, matrix 
quantum mechanics . However, we still hope that the limit c — > 1 can be found in the 
framework of GKM and, probably, by using of non-polynomial potentials. 

All these problems deserve further investigation and we hope to return to them else- 
where. 

We are grateful to L.Chekhov, A.Orlov, A.Sagnotti and A.Zabrodin for deep and very 
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useful discussions. 



Appendix A. Integrable hierarchies in the operator 
formalism 

Here we would like to describe some basic ingredients of the integrable hierarchies in the 
terms of the massless fermions. This approach was initiated in the series of papers [25] 
and it appears to be very fruitful for description of the general structure of the large 
variety of nonlinear integrable equations like KP, Toda etc. 
Let us define the fermionic operators on sphere 



i>(z) = £ ij k z k , r(z) = E r k z~\ (A.l) 
feez kez 

where fermionic modes satisfy the usual anti-commutation relations: 



iij k , r m } = hm , m = {r k , = 0. (a.2) 

The Dirac vacuum |0) is defined by the conditions: 

ip k \0) = , k < ; Vfe|0) = , k > . (A.3) 
We also need to introduce the "shifted" vacua constructed from |0) as follows 



\n) = < 



VVi-i • • • V'o |0), n>0 
y*...tl>*-i\0), n<0 
and satisfying the obvious conditions 



(A.4) 



ip k \n) = , k < n ; ip* k \n) = , k > n. (A. 5) 

From fermionic modes one can built the £/(l)-currents 



Jk = J2 Tp^i+k » J-k = Jk , k e Z + (A. 6) 



and define "Hamiltonians" 
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H(x) = ^2x k J k , H(x) = ^2y k J k , (A. 7) 

k=l k=l 

where {x k } and {y k } are half-infinite sets of independent time variables ("positive" and 
"negative" times correspondingly; in main text {x k } = {t k }, {y k } = {— t_ k }, k > 0) 
which generate the evolution of nonlinear system. 

Let g be an arbitrary element of the Clifford group which does not mixes the ip- and 
ip*- modes : 

9 =■ exp[^ A km ip k ^* m \ :, (A.8) 

where : : denotes the normal ordering with respect to the Dirac vacuum |0). Then it is 
well known [10,25] that 

r n (x,y) = (n\e H ^ge- R ^\n) (A.9) 

solves the two-dimensional Toda lattice hierarchy, i.e. is the solution to the whole set 
of the Hirota bilinear equations. Any particular solution depends only on the choice of 
the element g (or, equivalently it can be uniquely described by the matrix A km ). From 
the eqs.(A.2) one can conclude that any element in the form (A.8) rotates the fermionic 
modes as follows 

g^ k g~ l = ipjRjk , g^lg' 1 = r 3 R k }, (A.io) 

where the matrix R jk can be expressed through A jk (see [26]). We will see below that the 
general solution (A.9) can be expressed in the determinant form with explicit dependence 
of Rj k . In order to calculate r-function we need some more notations. Using commutation 
relations (A. 2) one can obtain the evolution of ip(z) and ip*(z) in times {x k }, {y k } in the 
form 

ip(z, x) = e H(x) 4>(z)e- H(x) = e^ z) i){z) , (A.ll) 

1>*(z, x) = e H{x) i>*{z)e- H{x) = e-^ z) ^*(z) ; (A.12) 
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if>(z, y) = e B{y) ^(z)e- R{y) = e f ^ '^(z) , (A. 13) 

i/>*(z,y) = e R{y) ^*{z)e- R{y) = e~^' z ~ 1 ^(z) , (A.14) 

where 



00 

£(x,z) = 5> fc * fc . (A. 15) 

fe=i 

Let us define the the Schur polynomials Pk(x): 



V\z\x\ = e^ z) = Pk{x)z k ; (A.16) 

k=0 

then from eqs. (A.11)-(A.14) one can easily obtain the evolution of the fermionic modes: 

00 

Mx) = e H ^ k e~ H ^ = £ i/>k- m P m (x), (A.17) 

m=0 

00 

Mx) ee e H ^r k e H{x) = E r k+m Pm(-x); (A.18) 

m=0 

00 

My) = e^Vfce"^ = £ ^k+mPmiy); (A.19) 



m=0 



^(y) = e^e"^ = £ (A.20) 

m=0 

It is useful to introduce the totally occupied state | — 00) (see definition (A. 4) in the limit 
n — > —00) which satisfies the requirements 

- 00) = , % e z . (A.21) 

Then any shifted vacuum can be generated from this state as follows: 

\n) = V'n-iV'n^-l - 00) . (A.22) 

Note that the action of an any element g of the Clifford group (and, as consequence, the 
action of e~ H<yV ^) on | — 00) is very simple: g\ — 00) ~ | — 00), so using (A. 18) and (A.19) 
one can obtain from eq.(A.9): 
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r n (x,y) = (-oo\...ip* n ^ 2 (-x)ip^ 1 (-x)gip n „ 1 (-y)^ n - 2 (-y)-\ - oo) ~ 

(A.23) 

~ det[(-oo|^*(-a;)^ i (-y)^" 1 | - oo)] 



\i,j<n-l 



Using (A. 10) it is easy to see that 



g*Pj(-y)g 1 = Y. p ™(-v)i , kRk, j+m (A.24) 

m,k 

and the "explicit" solution of the two-dimensional Toda lattice has the determinant rep- 
resentation: 

T n (x,y) ~ det H i+nd+n (x,y) <q , (A.25) 

where 

Hij(x, y) = Y. R k mPk-i{x)P m -j{-v) • ( A - 26 ) 

k,m 

The ordinary solutions to KP hierarchy [25] correspond to the case when the whole evolu- 
tion depends only of positive times {xk}; negative times {yu} serve as parameters which 
parameterize the family of points in Grassmannian and can be absorbed by re-definition 
of the matrix Rkm- Then r-function of (modified) KP hierarchy has the form 

r n (x) = (n\e Hix) g(y)\n) ~ det[^ J R fcj+ „(y)P fe _ i _„(x)] | iJ<0 , (A.27) 

k 

where g(y) = ge^ H ^ and 

R kj (y) = £ RkmPm-A-v) ■ ( A -28) 

m 

One can consider the reduction to the Toda chain hierarchy after imposing the condition 
on the element g [10] 

[Jk + Jk,g] = (A.29) 

which is equivalent to constraint 



[A + A -1 , R] = 
38 



(A.30) 



In this case 



gQ~ykJk — Q — ykJk^ — ykJkg gVkJk 

and r-function depends (up to the trivial factor) only on times {xk — yu}'- 

Tn {x,y) = e^ kXkVk (n\e H{x - y) g\n) . (A.31) 
The reduction (A. 30) has an important solution^ 

Rkm = Rk+m ■ (A. 32) 

In this case the matrix Hij defined by eq.(A.26) evidently satisfies the relations = Hi+j 
and 

(d Xk + d Vk )H i+ j = for any k < n — i , k < n — j (A. 33) 
due to the properties of Schur polynomials 

d Xn P k {x) = P k - n (x). (AM) 

The property (A. 33) certainly does not imply that the corresponding r-function depends 
only on difference of times because of restriction of values of k, but it restores correct 
dependence of times with taking into account of exponential in (A.31). 

Now let us establish the corespondence between matrices Rij (eq.(A.lO)), (eq.(A.26)) 
and Tij (eq.(2.12)), (eq.(2.2)). These equations transform into each other under charge 
conjugation of fermions ipk ^-k-i (^- e - vacua transform as follows: |n) —>■ \ — n)) and 
correspondence: = Rj_i ^ and = H^. 

19 Generally the solutions R n k — R n +k and R n k — R n -k are different, but for forced hierarchy, when 
r-function is the determinant of finite matrix, these two solutions are equivalent due to possibility to 
reflect matrix with respect to vertical axis without changing the determinant. 
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Appendix B. Fermionic representation of the matrix 
models 



Here we would like to discuss the fermionic language for (multi-) matrix models. Namely, 
we shall show that the r-function of the two-dimensional Toda lattice (A. 9) describes the 
whole variety of the (multi-) matrix models for some specific choice of the element g. 
Since one should reproduce the partition function of the matrix models from eq.(A.9), we 
shall deal with the forced hierarchies (see (1.5) and discussion in the sect. 2.1), i.e. 

Tn = , n < . (B.l) 
Therefore, it is reasonable to consider the point of the Grassmannian in the form 

9 = 9oP+, (B.2) 
where P + is the projector onto positive states: 

P + \n) = 6{n)\n). (B.3) 
There is exist a natural fermionic projector 



P+ =: exp[5>i#] : (B.4) 

i<0 

with the properties 



P + r- k = ^P-kP + = , k > ; (B.5) 
[P + M = [P + M=0, k>0; (B.6) 

Pi = P+- (B.7) 

The insertion of such projector into eq.(A.9) naturally leads us to conclusion that g 
should depends only on ipk and ipl with k > 0. We use the choice 



40 



5o =:exp|(^A(z,«;)V + (^;(^~ 1 )^^) -£^*j : 5 (B.8) 

where ip+(z) = J2k>oi J kZ k , — J2k>oi J k z ~ k an d 7 * s some contour of integration. In 

what follows we shall also use the projector 



P_=:exp[-5>i#]: (B.9) 

i>0 



with the properties 



P-^ k = r k P- = ,k>0; (B.10) 

[P_, V-fc] = [P_, W] = , k > ; (B.ll) 

P 2 = P_. (B.12) 

Now one should calculate the state 

(?oP + e^ fe) |n) • (B.13) 

It is easy to see that this state vanishes when n < 0. Indeed, using eqs.(A.4) and (A. 20) 
one can obtain that with n < the state 

e- s M\n) =p-n(-v) - Ci(-y)e-^|0) 

contains only negative modes ^*_ m (m > 0). Therefore the action of P + annihilates this 
state due to eq.(B.5). For n > using eqs.(A.4), (A. 19) and (B.6) we have 

P +e - R W\n) = iPn-ii-V) ... Vo(-y)P+e^ (y) |0). (B.14) 
Then we use the fact that 

P+e- B M\0) = |0) . (B.15) 

Proof of eq.(B.15). Let us denote 



41 



\y) = P + e- H M\0) . 

Then 

due to eqs.(A.3) and (B.5). Since \y) \ yk =o = |0), the eq.(B.15) is proved. 
Therefore, we have 

g P + e- S ^\n) = g j>(-y) . . . ^(-j/)|0) = 
1 r m 

= ^2~ ; J\A(z h w i )dz i dw i 'il) + (zi) . . .ip+(z m )x (B.16) 
ml J li=1 

xP_r + {w m v ) . ..r + {Oi> n -i{-y) ■ ■ -M-y)\Q) 

with using of the eqs.(B.8) and (B.9). Now we shall see that only the term with m = n 
gives a non-zero contribution in the infinite sum (B.16). Indeed, for m > n the state 
4'+{ w ^i) ■ ■ ■ 4>+( w i 1 )' l l'n-i{— y) ■ ■ • ^oC — 2/) |0) vanishes, because in this case some positive 
modes in ^(w^ 1 ) will reach the vacuum |0) and annihilate it. Vice versa, for m < n 
some positive modes in ipk{~V) wn l reach the projector P_ and due to eq.(B.lO) it is zero. 
Therefore, 



1 r n 

g P +e - H(y) \n} = — / T[A(z i ,w i )dz i dw i ip + (z 1 ) ...ijj + (z n )x 

xP_r + {w~ l ) . . . r+fa^n-ii-y) ■ ■ ■ M-V)\o)- 

Now we use the following proposition: 



r^n 1 ) ■ ..r + ( w ^n-i(-y) . ..m-v)\o) 



(B.17) 



;b.is) 



= AHex P [-^e(y,^)]|o) . 

Proof. Since the number of creation (w.r.t. to |0)) operators ipi(—y) equals to the number 
of annihilation operators ^(wj 1 ), then it is obvious that after normal re-ordering 

r+i^-r+^T^n-ii-y) ■ ..M-V)\o) = C( >nst ■ |o> 

and, consequently, 
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const = (oir+^n 1 ) ■ ■ ■ r+^T^n-ii-V) ■ ■ ■ M-y)\o) = 

= det[<0|^(«;r 1 )^-i(-i/)|0>] |y=i,...,„ 
and using eqs.(A.20), (A. 16) one can obtain 

(0|^(«;r 1 )^--i(-i/)|0>=^'- 1 e -€^ ; 

thus, 

n 

const = detK _1 e^ (2/w) ] = A(w) exp[- wj)}. 

3=1 

After substitution of eq.(B.18) into eq.(B.17) and using the obvious fact that P_|0) = |0) 
and ip_(zi)\0) = one can obtain 

<7 P+e-*M|n> = 

n (B.19) 

= -r / n^^<) e "^' Wi)d ^ A H^i)---^(^)l°)- 

Using one of the basic formulas [25] (which can be simply proved by bosonization tech- 
nique) 



rl>( Zl ) . . . ^(z n )\0) = A(z) exp[tf e(^)] |0>, (B.20) 

i=i 

where e(z;) is the vector with components €k(zi) = —z\ , we have the desired result: 

K 

g P + e- R M\n) = 

1 n n (B.21) 

= -f / n^(^ ) ^) e ^ (s/ * u,i)rf ^ A ( w ) A w ex p[^(E e (^)]i°)- 

n! - 7 Tj=l i=l 

Thus, finally we obtain: 

r n (x,y) = (n\e H ^g P + e- n ^\n) = 

1 , „ (B.22) 

- / A(w)A(z)l[A(z l ,w t )e^ x ' z ^ y ' m) dz i dw t . 



71! J 7 

Let us consider some particular cases. 



i=i 
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i) 7 — ( — °°> +°°) , A(z,w) = S(z — w). Then one can recover the case of Hermitean 
one-matrix model. 

ii) 7 is the small circle around the origin in the complex z-plane, and A(z,w) = 
5{z — w^ 1 ). Then g — 1 and r n (x — y) is the r-function for the symmetric unitary 



2m z 
model. 

hi) 7 = (— oo, +oo) , A(z,w) = e zw . This is Hermitean two-matrix model. 

iv) 7 = (— oo, +oo) . Let us denote Zi = zf^ , Wi = zf 1 , x = —t^ , y = and 

A{zf\ = - / II ex P{" E + E (B-23) 

Jo ° j=l 3=2 3=1 

Then we obtain the partition function for Hermitean p-matrix model which originates 
from the matrix integrals 

f DZ^K..DZ^exp{-TrY,^\Z^)+TrJ2Z^Z^ +1) }. (B.24) 

3=1 3=1 

The only point is that in this case g is parameterized by the set of additional times 

j — 2, . . . ,p — 1. One can obtain the multi-matrix model with time-independent g in the 

context of the multi-component Toda lattice hierarchy. 

Determinant representation. 

Using eqs.(A.4), (B.7) and the fact that [P + ,g ] = 0, we have: 

r n (x,y) = (0\r o ...r n -ie H{x) 9oP + e- n ^^ n . l ..^ \0) = 

= (o\e H ^r (-x) . ..r n -i(-x)P + 9 P + ^n-i(-y) ■ ..M-y)e- R{y) \o) ■ 

Since ip*(—x) and ipi(-y) contain only positive modes (see eqs.(A.18) and (A. 19)), due 
to (B.6) and (B.15) one can obtain 

r n (x,y) = (o\i/>*(-x) • • -i>n-i(- x )9oi>n-i{-y) • ••-0o(-y)|o) = 

(B.25) 

det[(0\^(-x)g o ^(-y)\0)} \i, j=0 ,...,n-i ■ 
The same arguments used for transition from eq.(B.16) to eq.(B.17) when applied to 
eq.(B.25) lead to conclusion that only linear term in A(z, w) contributes, so using eq.(B.lO) 
we have 
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(OK(-x)^-(-y)|0> = / A(z,w)dzdw(0\^(-x)M^)P~r + ^~ 1 )M~yM = 
= f z t w 3 A{z,w)e^ x ' z) -^ y ' w) dzdw = di{-d y ) j f A{z,w)e^ x ' z ^ (y ' w ^dzdw 

(B.26) 

and, finally, one can obtain the expression for r-function in the determinant form: 

r n (x,y) = det^-dy,)* f A(z,w)e^~^dzdw] | M=0 ,...,n-i . (B.27) 

J-y 

Again, the consideration of particular choices of A(z,w) (see discussion below eq.(B.22)) 
leads to representation of r-function for Hermitean, unitary etc. models in the determi- 
nant form. 
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